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Abstract 



We apply the multicanonical technique to the three dimensional 
\ dynamical triangulation model, which is known to exhibit a first order 

phase transition with the Einstein-Hilbert action. We first clarify the 
first order nature of the phase transition with the Einstein-Hilbert 
action in several ways including a high precision finite size scaling 
analysis. We then add a new local term to the action and confirm the 
conjecture made through the MCRG technique that the line of the 
first order phase transition extends to the expanded phase diagram, 
ending at a point. Fractal dimension at the end point is measured to 
be around three up to the present size. 
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1 Introduction 



Quantum gravity is one of the most important topics in particle physics. The 
attempts to construct it within ordinary field theories confront a difficulty 
due to the fact that the theory is perturbatively unrenormalizable. This does 
not mean immediately that the attempts fail, since it might be possible to 
construct it nonperturbatively as is the case with the 3D nonlinear sigma 
model. The lattice regularization provides a natural approach in this con- 
text. In the case of quantum gravity, since the dynamical variable is given 
by the metric of the space time, how to discretize the space time is itself 
a nontrivial problem. Above all, it seems difficult to preserve the general 
coordinate invariance manifestly. The success of the dynamical triangulation 
model in two dimensions []T| summarized by its equivalence to the Liouville 
theory [0] in the continuum limit shows that the model provides a natural 
framework to study quantum gravity nonperturbatively within ordinary field 
theories. Although the restoration of the general coordinate invariance in the 
continuum limit is only supported by the success in two dimensions, we might 
consider that the summation over all the triangulation with the equal weight 
forms the basis of the restoration. The situation with the other lattice reg- 
ularization of quantum gravity known as the Regge calculus is very obscure 
i- 

The above point of view motivated intensive studies of the dynamical 
triangulation model in higher dimensions through Monte Carlo simulation. 
It turned out, however, that with an action obtained by naive discretization 
of the Einstein-Hilbert action, the phase transition observed is of first order 
both in three dimensions and in four dimensions ||5|, |^. The signal of 
the first order phase transition is much weaker in 4D than in 3D. We note 
that this does not necessarily mean that it is easier to construct a continuum 
theory out of this model in 4D than in 3D, and we consider that this is merely 
related to the existence of singular vertices in the crumpled phase in 4D J^, § . 
By adding a local term which suppresses the singular vertices, we observe a 
first order phase transition in 4D, which is even stronger than in 3D. Indeed, 
from the viewpoint of perturbative unrenormalizability, the situation in 4D 
must be severer than in 3D. The e-expansion of (2 + e)-dimensional quantum 
gravity [Q suggests the existence of a non-trivial fixed point, and this result is 
of course more reliable for smaller e naively. We therefore concentrate on 3D 
dynamical triangulation model in this paper. Three-dimensional quantum 
gravity is worth being studied in itself, since it can be also regarded as the 
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first quantization of membranes, which is considered to be a fundamental 
object in the M-theory [|1^ . 



In lattice theories in general, it is often preferable {e.g. improved actions 
to reduce the finite lattice spacing effects) or even necessary {e.g. the Wil- 
son term to eliminate fermion doublers) to modify the lattice actions from 
that obtained by naive discretization of the action in the continuum. We 
may therefore hope to obtain a second order phase transition in the present 
case by adding a local term to the naive Einstein-Hilbert action. Once the 
continuum limit can be taken, the continuum theory constructed is expected 
to be independent of the details of the constructions thanks to universality. 
Adding some matter fields or gauge fields is another possibility to obtain 
a continuum theory, though it changes the theory to be constructed. We 
stick to pure gravity as the simplest case in this paper. 

We perform a numerical simulation of this model using the multicanonical 
technique [0, which is powerful especially in the presence of large poten- 



tial barriers between some local minima, as in the present case in which 
the system undergoes a first order phase transition. We measure the quan- 
tity which corresponds to the 'temperature' in the microcanonical ensemble, 
which enables us to extract various important thermodynamical properties 
of the model very clearly. We first clarify the first order nature of the phase 
transition with the Einstein-Hilbert action in several ways including a high 
precision finite size scaling analysis. We next add a new term to the action, 
which is motivated by a possible contribution from the measure for the path 
integral over the metric. The modified action has been first studied in four 
dimensions fl^. Recently it has been studied in three dimensions by the 



Monte Carlo renormalization group and a possibility was suggested that the 
first order phase transition ceases to exist at a finite coefficient of the new 
term [Q. We confirm this possibility by investigating the expanded phase 
diagram []. 

The phase transition is expected to become second order at the end point 
of the phase transition line. We therefore attempt to observe the continuum 
physics at this end point. We measure the boundary area distribution, whose 



counterpart in 2D is known to show a clear scaling behavior [16, 1/]. We 
indeed see a reasonable scaling behavior. 

The paper is organized as follows. In Section S, we explain the mul- 



^ While this paper was being completed, a preprint appeared |15[ , in which the authors 
address the same issue but with a conclusion which is different from ours. 
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ticanoical techniques and introduce the quantity which corresponds to the 
'temperature' in the microcanonical ensemble. In Section ^ we show the 
numerical results with the Einstein-Hilbert action, and clarify the first order 
nature of the phase transition. Above all, we observe a clear finite size scal- 
ing, which is consistent with the first order phase transition. In Section ^, 
we introduce the model with the new term. In Section ^ we show the results 
with the modified action and confirm that the first order phase transition line 
extends to the expanded phase diagram, ending at a point. In Section ||, we 
show the results for the boundary area distribution at the end point, which is 
consistent with fractal dimension three. Section ^ is devoted to conclusions 
and discussions. 



2 Multicanonical technique and "microcanon- 
ical inverse temperature" 

The naive discretization of the Einstein-Hilbert action in the dynamical tri- 
angulation model in three dimensions gives Q 

Seh = -hqNo + K3N3 , (1) 

where A^^o and N3 are the number of vertices and 3-simplices respectively in 
the simplicial manifold. Kq corresponds to the inverse of the gravitational 
constant, while k,^ corresponds to the cosmological constant. Since we con- 
sider ensembles with fixed N^, we use only the first term So = —kqNo in the 
following. 

The partition function of the system with fixed A^s is given as follows. 
Z{ko;Ns) = J2 exp(-S'o) 

{T{N3)} 

= ^P(Aro, /to; ATg), (2) 

P{No,^o;Ns) = n{No;Ns)exp{KoNo). (3) 

The summation in the first line of (0) is taken over all possible triangulations 
of 3-dimensional simplicial manifolds with sphere topology that contain A3 
3-simplices. n{No] N3) is the number of the possible triangulations with a 
further restriction that they contain Aq vertices. P{Nq, hq; N-^) represents 
the node number distribution for a canonical ensemble with fixed kq. In the 
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following we suppress the argument in the P{Nq, Kq; N^) and u^Nq; N^,), 
unless necessary. 

As in all of the recent numerical works on higher dimensional dynamical 
triangulation model, we use the so-called {p,q) moves, which are proven to 



be ergodic [0, to simulate this system. Actually, is not kept fixed in 
the moves. We therefore have to allow a small fluctuation of around a 
target with an appropriate Gaussian potential for N^. Measurement has 
been done only when the agrees with the target N3. We refer to it as a 
simulation for fixed A^3 for simplicity. 

The system is known to undergo a first order phase transition at some 
critical kq. Near the critical point, P{Nq, kq) exhibits a double-peak struc- 
ture. As one increases the system size N^, the valley between the double 
peaks becomes deeper and deeper, which makes the ordinary canonical sim- 
ulation increasingly difficult, due to the long Monte Carlo time required to 
travel back and forth between the two peaks for sufficiently many times. 
We perform the multicanonical simulation to overcome this problem. For 
technical details, we refer the readers to Ref. |T^. The simulation is per- 
formed by replacing the Boltzmann weight exp(/toAi"o) with some function of 
No, which is determined iteratively so that the distribution of Nq becomes 
as constant as possible within some window of No. If the window is wide 
enough to cover the whole support of the distribution function of A^o for a 
canonical ensemble with some kq, one can reproduce the expectation values 
of any observables in the canonical ensemble by taking account of the effect 
of replacing the Boltzmann weight in updating the system. The idea is essen- 
tially a generalization of the histogram method. If one extends the window 
by fine-tuning the Boltzmann weight, the simulation includes information of 
canonical ensembles with various kq in some finite region. This enables, for 
example, a precise identification of the peak of the node number suscepti- 
bility as a function of kq, which is essential in the high precision finite size 
scaling analysis. 

Another advantage of the multicanonical simulation is that it allows us 
to extract the entropy as a function of A^o defined by 



SiNo) = log niN,) , (4) 



or its first derivative 



KiNo) = -{S{No + 1) - SiNo)}. (5) 
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K{Nq) is the analogue of the inverse temperature (3 = 1/T defined in the 
microcanonical ensemble if we regard —Nq as the energy. We therefore refer 
to K{No) as the microcanonical inverse temperature. Strictly speaking, the 
simulation allows us to extract the entropy S{No) only up to an additional 
constant independent of No. Once S{Nq) is known, one can calculate the 
node number distribution P{Nq, Kq) in the canonical ensemble for arbitrary 
through 

P{No, Ko) = exp{ft:oiVo + S{No)}. (6) 

From this, we have 

logP(iVo + 1, Ko) - logP(iVo, Ko) = /to - K{No), (7) 

which means that if kq — K{No) is positive (negative) for some No, the node 
number distribution P{No, kq) increases (decreases) as A^o increases. This 
is illustrated in Fig. |I]. The upper part shows the microcanonical inverse 
temperature K{Nq) as a function of Nq, with three horizontal lines K = kq 
corresponding to kq = Kq, k"^ and Kq. The lower part shows the node number 
distribution P{No, ko) for each of the three values of kq- The intersections of 
the curve of the microcanonical inverse temperature and the horizontal lines 
give the extrema of the node number distribution for each kq. 

When Kq = K^, there are three intersection points: = K{Nq'^), 
K{N^^^), and K{N^^^), where N^'^^ < N^^^ < N^^K The node number 
distribution has the double peaks at Nq'^ and A^^q^-*. Nq'^ gives the place 
of the local minimum. From (0), the ratios P{N^^\ k'^)/P{N^^\ k"^) and 
P{Nq , K^)/ P[Nq , K^) are given by the exponential of the area of the light 
and heavy shaded regions in Fig. respectively. Since we have chosen 
so that the area of the two regions is equal, the peaks of the node number 
distribution have the same height. 

In general, when the microcanonical inverse temperature has a local min- 
imum and a local maximum, the node number distribution shows a double 
peak structure for the canonical ensemble with set to the value between 
the local minimum and the local maximum. The quantity, therefore, serves 
as a useful probe to clarify the fate of the first order phase transition when 
one modifies the action from the Einstein-Hilbert action. 



6 



Microcanonical inverse temperature 
and the node number distribution 




3200 3400 3600 3800 4000 4200 



No 

Figure 1: The upper part shows the microcanonical inverse temperature 
K{No) (heavy sohd curve) with three horizontal lines K = kq corresponding 
to Ko — Kq, kIq and Kq. The lower part shows the node number distribution 
P{No,K,o) for each of the three values of kq. 
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microcanonical inverse temperature 

Kws Ng/N, 




3.6 I ' ' ' ' ' ' 1 

0.25 0.27 0.29 0.31 

Figure 2: The microcanonical inverse temperature K{No] N^) is plotted 
against the node number A^o normalized by the system size for = 3, 500 
to 40,000. 

3 The results for the Einstein-Hilbert action 



In this section, we show the results for the Einstein-Hilbert action and clarify 
the first order nature of the phase transition by the measurement of various 
thermodynamical quantities. 

We first show the microcanonical inverse temperature K{No; N^) for vari- 
ous system size from = 3, 500 to 40, 000 in Fig. |[ Since the microcanoni- 
cal inverse temperature has a local minimum and a local maximum, the node 
number distribution shows a double peak structure at some region of kq. We 
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choose the kq with which the height of the peaks is the same. We denote such 
'^o by Kq^. The determination of Kq'^ can be done easily by referring to the 
microcanonical inverse temperature as is explained in the previous section. 
Fig. ^ shows the results for the node number distribution at k'q' with the 
system size from N3 = 3, 500 to 40, 000. The logarithm of the ratio of the 
height of the peaks to that of the valley gives the interfacial energy. (See 
Ref. \T^\ for further analysis on this quantity.) Note that the ratio is about 
10^^ for N3 = 40, 000. One can easily understand that ordinary canonical 
simulations cannot be used in this situation. 

In Fig. ^, we plot the mean node number in the canonical ensemble 
with against kq- We replot the microcanonical inverse temperature with 
kq = K{Nq] N^) for comparison. The mean node number shows a sharp 
increase in the critical region of hq for a sufficiently large system size, which 
is another way of observing the signal of the first order phase transition. 
The data curve of the mean node number coincides with the data curve of 
the microcanonical inverse temperature except in the critical region, as is 
expected. 
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yV/A/g vs Kg and K 




Figure 4: The mean node number normalized by the system size N3 is plot- 
ted against k,q by the solid line. The microcanonical inverse temperature is 
replotted by the dashed line for comparison. 
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Figure 5: The node number susceptibility xii^o'i ^s) is plotted against kq for 
the system size from = 3, 500 to 40, 000. 



In Fig. ^ we show the result for the node number susceptibility xi^o] ^s) = 
{{N^) - {Noy)/N3. We call the kq which gives the peak of the susceptibility 
for each N^, as the pseudo-critical point KqIN^). In ordinary statistical sys- 
tems, the peak height of the susceptibility grows as a power of the system 
size. If the peak is due to a first order phase transition, the internal energy 
distribution shows a double peak structure as the one in Fig. |, and the 
tunneling between the two peaks causes a linear growth of the peak height 
of the susceptibility. In Fig. | we plot the peak height of the susceptibility 
against the system size. The data for the four large N^^s can be fitted very 
nicely to the formula 

xiK^oiNs); Ns) = pi + P2N3 , 

with pi = —1.8(2) and p2 = 3.20(7) x 10^'^, which is consistent with the finite 
size scaling for the first order phase transition. 

The pseudo-critical coupling constant KqIN^) shifts to the right with no 
tendency of convergence as one increases the system size up to A^3 = 40, 000. 
This behavior should be contrasted with ordinary statistical systems, for 
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peak height of susceptibility vs 

15.0 




Figure 6: The peak height of the susceptibihty is plotted against A^3. The 

straight Une show the hnear fit x/^s = Pi + P2N3, which is consistent with 
the first order phase transition. The of the fit is /dof — 0.3. We use 
only the data for the four large A^s's. 
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which the dependence of the pseudo-critical couphng on the system size V 
obeys the scahng relation such as 



(8) 



The unusual shift of the pseudo-critical coupling constant has been discussed 
earlier in the four- dimensional dynamical triangulation model pO|. 



4 The modified action 

In this section, we modify the Einstein-Hilbert action and attempt a 
search for a second order phase transition point in an expanded phase dia- 
gram. We consider the following action. 

S = 5o-/iM, (9) 
M = ^log[o(t;)/4] , (10) 

V 

where is the new coupling constant 0, and o{v) is the number of 3— simplices 
that contain the node v. The summation in (0) is taken over all the nodes 
in the manifold. This term |T^ can be considered as a possible contribution 



from the path integral measure II^ g{x)^^'^, since o{v) can be interpreted as 



the "local volume", which corresponds to y g{x) in the continuum. 

In the following we consider the canonical ensemble with fixed given 
by the partition function : 

Z('"°")('^o;/U;iV3) = Yl exp{KoNo + fiM), (11) 

mNs)} 

= ^;2(iVo;/i;iV3) exp(KoiVo), (12) 

where z{Nq] /x; N^) is the partition function for fixed A^o; which can be given 
by 

z{No;fi;Ns)= ^ exp(/iM) , (13) 

{T{N3,No)} 



^Note that the u defined here is the same as the one in Refs. 0, The definition 
of fi used in Rcf. pi[ corresponds to — /i in tliis paper. 
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where the summation is now taken over all possible triangulations that con- 
tain N3 3-simplices and Nq vertices. As before, we suppress the argument 
N3 in z{No] fi] iVs) unless necessary. In this paper, we concentrate on the 
negative fi region, in which the first order phase transition line is predicted 
to end at a point through the MCRG analysis in Ref. [ff4 



5 End point of the first order phase transi- 
tion 

Let us define the (generalized) microcanonical entropy 

S{No;fi) = log z{No;fi) (14) 

and its first derivative 

K{No; fi) = -{S{No + I; fi) - S{No; /i)}, (15) 

which we keep on referring to as the microcanonical inverse temperature. 

In Fig. 1^, we show the microcanonical inverse temperature K{No] n) for 
several negative values of /i. As was explained in Section H, the double peak 
structure in the node number distribution, which is a signal of the first order 
phase transition, can be probed directly by looking at the microcanonical 
inverse temperature K{Nq] /i). 

One can see that the local minimum and the local maximum observed 
for /i = merge at /i* = —0.65, beyond which the K{No; /x) becomes a 
monotonous function of Nq. This suggests that the first order phase transi- 
tion ceases to exist at /i*, beyond which it becomes a crossover. 

For each /i > /i*, we define the pseudo-critical coupling Kq as the position 
of the peak of the node number susceptibility Xo{f^o, lA- The end point of the 
phase transition line is given by /i* = —0.65, Kq = 5.00 for = 10, 000, and 
/i* = -0.80, Kg = 5.515 for A^'s = 20, 000. In Fig. |, we plot the free energy 
/(A'o; yu; /^o) = —S{No; fi) — kqNq as a function of Nq near for /i = and 
H = fi*. The two local minima observed for fi = merge into one at yU = /i*. 
The free energy curve at {kq, fi*) has a wide flat bottom, which suggests the 
existence of a massless mode. 
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with |i term 



^ 5.6 

51 



fc< 5.3 



5.0 





^=-0.50 


1 


r 




^=-0.58 








^=-0.62 




1^7 - 




^i*=-0.65 








^=-0.70 


rr 






H=-0.80 








^=-0.85 








^=-0.90 


„/■■'•' r 






^=-1.00 ^ 


^ /••.-■'"' -'" 'ZiT' 

,-i"rh 

















— — — """^r 





0.29 0.30 0.31 0.32 

Figure 7: The microcanonical inverse temperature K{No; /x) as a function of 
A^o for several /.t's with the system size A^3 = 10, 000. 
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Free energy with \i term 
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Figure 8: The free energy /(A^o; yU, kq) as a function of Nq near for yU = 
and fi = fi* = —0.80 with the system size = 20, 000. 



6 Fractal structure at the end point 

In order to examine a possible continuum hmit at the end point, we measure 
the boundary area distribution p{A,r), which is the number of boundaries 
with the area A at the geodesic distance r. The corresponding quantity in 
two dimensions called the loop length distribution is calculated analytically 
and is found to possess a continuum limit \W\. The scaling behavior expected 



from this result has been correctly reproduced by a numerical simulation |T7 
In Fig. we plot the boundary area distribution at the end point of the 

first order phase transition line as a function of x = Ajr^. One can see a 

reasonable scaling behavior. The power of r in the scaling variable x = Ajr^ 
implies that the fractal dimension is 3. We would like to remark that 

the scaling behavior we observe is much better than the one that has been 



claimed to exist with the unmodified action [22 
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Boundary area distribution 

N,=20,000, \i=-0.80 K„=5.515 

10^ I ^ 




x=AJr 



Figure 9: Boundary area distribution p{A^ r) at the end point of the first 
order phase transition fine with the system size A^a = 20, 000. 
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7 Conclusions and discussions 



We applied the multicanonical technique to the three dimensional dynamical 
triangulation model. The microcanonical inverse temperature served as a 
useful tool to clarify the first order nature of the phase transition observed 
for the Einstein-Hilbert action and to investigate its fate in the expanded 
phase diagram with a new coupling constant. 

Our results about the modified action suggest that the first order phase 
transition becomes second order at the end point, where one can construct a 
continuum theory. A natural question to be asked is what the continuum the- 
ory can be. One may suspect that we cannot construct any continuum theory 
in three-dimensional quantum gravity, since we have no physical degrees of 
freedom. This argument is too naive, however. An example of field theory 
which has no physical degrees of freedom that is still completely well-defined 
is two-dimensional Yang-Mills theory. 

The fact that we have to fine-tune two parameters to obtain the con- 
tinuum limit implies that there are two relevant operators around the fixed 
point. A possible interpretation of the continuum theory is therefore the 

gravity. The observation that the fractal dimension at the fixed point is 
approximately three is also suggestive of this interpretation. The fractal di- 
mension we extracted, however, is still preliminary, and the real value, which 
could be obtained by increasing the system size might well be larger. 

In Ref. |T^, the authors pointed out that the system in the negative 



/i region might be suffered with large finite size effects. Taking this into 
account, they claimed that the line of first order phase transition continues 
to ^ = fi'^ and kq = oo and that the system is in the crumpled phase for 
all Kq when /x < /i"^. Needless to say, simulations with larger system size are 
necessary to discriminate between these possibilities. 

In Ref. 0, we examined a different modified action in four- dimensional 
dynamical triangulation model. The added term was 

S^ = uJ2[o{v)-5]\ (16) 

V 

We saw that the phase transition turned into a cross over even for a small 
positive u and the system was observed to be always in the branched polymer 
phase. Let us write the term to be added as J2v fio{v)). Note that a constant 
term and a linear term in f{x) can be absorbed by the redefinition of kq and 
K3. The quadratic term in (|T6[) might be too strong since it dominates the 
behavior of f{x) in the large x region. 
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In Ref. |T^, another possibility has been suggested that the hne of the 
first order phase transition ends at a point with positive /i. It is of interest 
to clarify this possibility as well using the technique we used, although the 



simulation must be much harder due to the small acceptance rate [0. 
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